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7^ , Abstract: A handy and extensively used kinematic variable in LHC analyses, especially for 

golden SUSY signals of multijets plus missing energy, is the Effective Mass, Meg = ^ Ipt I + br "'^'^ I ■ 
^ Empirically, the value of A/off at which the histogram of events has a maximum is correlated with 

•/"J the SUSY spectrum, Afoff jj^j^^^ ~ 80% Msusy, where Afgusy is essentially the sum of the masses of the 

SUSY particles initially created. In this paper we explain the reason for such strong correlation, 
pointing out the cases where the correlation is not good. Besides, we propose a new variable, 
the Effective Transverse Energy, which shows an even better and more direct correlation 

^T^lmax ~ -^susy, ^nd is independent of the procedure followed to identify the jets. and Afcfi 

T-H are complementary variables, rather than competitors; and plotting histograms in both can be useful 

J> to cross-check the results, allowing a better and more robust identification of M^nsY- The extension 

of this procedure to other scenarios of new physics (not necessarily SUSY) is straightforward. 

Keywords: Beyond Standard Model, Supersymmetry Phenomenology, Supersymmetry 
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1. Introduction 

The LHC is already probing new physics beyond the reach of past experiments. There are 
two main questions to address: 1) Is there any signal of New Physics (NP)? and 2) In the 
positive case, which NP is it? In order to optimize the answer to these questions there is an 
intense activity to explore assorted strategies for the search of NP. The task is challenging, 
due in part to the fact that LHC data, though very rich, are not as clean as those from an 
e"*" e~ collider. 

Supersymmetry (SUSY) is one of the few candidates for physics beyond the SM that is 
really well motivated from the theoretical point of view and allows to perform detailed cal- 
culations and thus realize precise predictions. Indeed SUSY has been the most extensively 
studied candidate for new physics in the last decades and the first LHC analyses are using 
SUSY as a paradigmatic scenario of new physics to present their constraints on physics 
beyond the SM. Of course, this does not mean that SUSY is really there, but clearly is a 
most serious scenario to be considered in the light of the LHC. 

In the framework of SUSY and under the assumption that superparticles are light 
enough to be produced and detected at LHC, one can expect a typical signal of multijets 



- 1 - 



and large missing energy with or without leptons. Searches using channels with leptons 
in the final state allow quite clean reconstructions. Looking for end-points and using in- 
genious kinematic variables [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11], one can determine in some 
cases the mass of the sparticle produced in the decay chain. On the other hand, the chan- 
nel with zero leptons tends to be the golden one producing signals beyond the Standard 
Model. However, here the analysis becomes much more complicated and the correct iden- 
tification and reconstruction of the properties of the final states is quite tough. Besides, 
the separation of the decay chains is normally not possible. 

In this context, the use of appropriate kinematic variables and the choice of optimal 
cuts are instrumental to maximize the SUSY signal from the Standard Model background. 
But, for the multijet channel, finding the optimal choice is a difficult task. In ATLAS 
analyses the Effective Mass variable (Mgff) plays a very important role [12, 13, 14]. It is 
defined as, 

4 

M!:S = Y.w\ + \prn, (1-1) 

where prp is the transverse component of the momentum of the four hardest jets [thus the 
superscript (4)] and p^^"^*^ is the missing transverse momentum. 

It was found in [15] (from now on HPSSY) that for the Constrained Minimal Su- 
persymmetric Standard Model (CMSSM) there is an outstanding correlation between the 
value of at which the histogram of events has a maximum and the supersymmetric 
mass Msusy defined as the minimum of the squark and gluino masses, 

^fsf ^""^ = min{m,-, m^}. (1.2) 

f4) fHPSSY) 

More precisely, HPSSY showed that systematically Mj^ [max — 1.9 Msusy for a 
large collection of CMSSM models with assorted values of the initial supersymmetric pa- 
rameters. This is very remarkable. However there are some caveats: 



1. As admitted by HPSSY the definition of the supersymmetric mass as in eq.(1.2) 
was rather arbitrary. Although typically the production of a pair of the lightest 
supersymmetric particles (among squarks and gluinos) is favoured, the production of 
other pairs, like squark- gluino, can be very frequent too. A refined definition of Msusy 
was proposed in [16], namely the sum of all supersymmetric masses weighted by the 
production cross section (normalized to one) of each one. In this way, the correlation 
behaves even better, though since the new Msusy is slightly larger the above number 
1.9 decreases somewhat. 

2. On the other hand, M^^^ is also smaller than the total effective mass, M^g . A 
more compelling definition, from the theoretical perspective, would be to sum in 
eq.(l.l) over all jets (surviving certain cuts); though this might be problematic from 
a practical point of view. 
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3. So far there is no theoretical explanation for the strong correlation found between the 
effective mass and the supersymmetric mass. In ref.[15] it is presented as an empirical 
fact. Certainly, for events with large transverse momenta, one would expect the 
effective mass to be of the order of sum of the masses of the supersymmetric particles 
initially created. However, it is not clear why the maximum number of events is 
always reached at a value of the effective mass so strongly correlated (but slightly 
smaller) to that sum. 

4. As a matter of fact, it was shown in [16] that for more generic MSSM models (beyond 
the CMSSM) such correlation fails in many instances. Again, it would be very useful 
to know the reason for that, in order to avoid flawed analyses and improve general 
strategies. 

In the present paper we will offer an explanation for the above points (3) and (4) . This 
understanding will allow us to propose a new kinematic variable, alternative to the effective 
mass, which shows an even better correlation with the supersymmetric masses, and can be 
used as an alternative or complementary check in multijet studies at LHC. The discussion 
and proposal presented here can be easily extrapolated to other scenarios of New Physics 
different from SUSY. 

In section 2 we introduce some relevant kinematic concepts and establish an analogy 
between the production of the W boson and that of Supersymmetric particles in hadron 
colliders. In section 3 we explain the empirical strong correlation between M^q and Mgusy 
In section 4 we propose a new kinematic variable, which shows an even more robust 
correlation with Mgusy. In section 5 we test the efficiency of M^s and and compare 
their behaviour using signal simulations at the LHC. 

2. Strategy for the kinematic analysis 
2.1 W boson production in pp colliders 

We start by briefly reviewing the W boson production and decay in a pp collider, since 
it is the simplest case involving visible and invisible particles, and it is therefore a useful 
guide to introduce some notation and some relevant kinematic concepts which will be used 
later on. 

In a pp collider the main process for W boson production is through quark- antiquark 
annihilation. Of course, the laboratory (LAB) system of reference does not coincide in 
general with the center-of-mass (CM) one, since it is affected by boosts, mainly along 
the collision line (z). In addition, there can be less important boosts along a transverse 
direction due to other effects; in particular, the quarks in initial state may radiate soft 
gluons (some of them energetic enough to be detected as jets), but for the moment we will 
ignore these effects. 

The W boson mass is measured by analyzing its leptonic decay — )• — )• ez^, see 
Fig.l. From the energy and momentum of the electron and the (anti)neutrino, one can 
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Figure 1: Feynman diagram of W-boson production and decay at hadron colliders. 

derive the invariant mass of the W boson, 

M^ = Ml, = {E'^ + E'^f-ip'+p'^f (2.1) 
= ml + ml + 2[E^E^cos\i{y^ -y") -Pt'Pt] , 

where is the transverse component of the 3-momentum of the a— particle, while Ef^ and 
are respectively its transverse energy and rapidity, defined as 

{E^f =ml + ip^j.) ' = {E'^f - [plf , y^ = \ In (|^) • (2.2) 

Note that Eqp is invariant under boosts in the z— direction. E"" and can be written in 
terms of E^ and as 

E'' = E^ cosh y'^ , p^ = E^ sinh y"". (2.3) 

Since it is not possible to measure the neutrino momentum, one cannot directly deter- 
mine the W invariant-mass from the experiment using eq.(2.1). Instead, one can consider 
the differential cross section of production of an electron-neutrino pair. 

Let us start working in the CM system. Calling 6, (p the polar and azimuthal angles of 
the electron 3-momentum, the corresponding differential phase-space factor is proportional 
to the solid angle element da oc d?il = d cos Odcp, where the integration over cp can be 
performed in a trivial way using symmetry around the z— axis. Now, one can consider the 
distribution of cross section for different values of p^, i.e. 

. dcos9 . „ ,^ ^, 

^ Hp^^^ ' ^ 

where the Jacobian factor reads 

dcos^ —p^\pe\ 



-2E pf^\pe 



^£;2 _ (^e + El^y^/E^ - {E^ - E-y ' 



(2.5) 
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with E = E'^ + E'^ denoting the total energy at CM {E = M\y if the W is produced on- 
sheU). Of course the neutrino and electron masses are completely negligible in this context, 
but for future convenience we will maintain them in the expressions. 

The phase-space factor (2.5) shows a pole at E^ + E!^ = E. Denoting 

£t = E^ + E!^ = ^(pe,)2 + ml + ^{p-^y + mi (2.6) 

(which is ~ 2p^ if nie and rriu are neglected), one can change variables — t- £t, introducing 
a trivial Jacobian factor. So in the £t variable the pole occurs at 

frlpoie = Ec:,Mw . (2.7) 

In consequence, the histogram of events in the £t variable should show a peak at £t = E ^ 
Mw followed by an abrupt fall to zero, since £t ^ E, . This behavior is slightly softened 
by the fact that the W has a non-vanishing width. 

The previous analysis was done at the CM system, but its extension to the LAB one 
is straightforward. First of all, since £t is invariant under z— boosts, the presence of the 
pole in at the CM energy (~ Mw) holds in the LAB system assuming there is no 
net transverse momentum, which is a sensible approximation. Actually, this is a successful 
strategy to identify M^y in this context. Of course, in practice there can be a net transverse 
momentum if the collision is not completely along the z— line, due mainly to initial state 
radiation. To incorporate this effect to the analysis note that, in CM frame, the electron 
and neutrino momenta for the events at the pole have no longitudinal components, so that 
Ej, = E^, E^ = E'^, and £t = E^ + El^ coincides with the total energy. These features 
obviously hold after a boost in a transverse direction. Consequently, in the boosted (LAB) 
system the pole in the cross section still occurs at 

^Tlpoie = E = ^Elj^, + {p^+p-f , (2.8) 

where Ecm = Mw. (Of course the previous expression is invariant under further z— boosts.) 
If desired, the effect of a net transverse momentum of the W can be extracted by defining 
a "transverse mass", Mt, as 

M| = 4 - {pt' + p^f = {E^ + E!^f - {pt" + Pt? 

^2|p/||p/|[l-cos(0^-0'^)], (2.9) 

where in the last expression we have neglected the electron and neutrino masses, which is 
allowed in this context. Note that Mt is invariant under any boost and satisfies Mt < 
Afinv, so the cross section does present a pole at Mt = Mw- The complication is that 
the determination of p^ + pjf is not as clean as one would like (although normally is a 
subdominant term). Under the assumption of a perfect longitudinal collision, this term 
vanishes and Mt coincides with £t- 
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Figure 2: Diagram of production and decay of supersymmetric particles at LHC. 
2.2 Pair production of SUSY particles at LHC 

Now we consider our scenario of production of supersymmetric particles (in pairs) at the 
LHC. In most cases the supersymmetric particles produced are squarks and/or gluinos, 
which decay along jets through diverse channels, plus one lightest supersymmetric particle 
(LSP, usually a neutralino) for each SUSY particle. Therefore the final state mainly consists 
of jets (plus possibly some leptons) and two LSPs (plus possibly some neutrinos). This is 
schematically represented in Fig. 2, where we have not shown the fact that the decays ocurr 
along two chains, one for each initial supersymmetric particle. Normally the SUSY particles 
are created close to at-rest in the CM system, since producing them with substantial 
momentum amounts a high price from the parton distribution functions. Along the paper 
we will often use this approximation. 

To connect with the strategy of the previous subsection, we can imagine that this 
process consists of the production of two pseudo-particles, J and X, which contain the 
jets and the two invisible LSPs respectively, as shown in Fig. 2. Hence the momentum 
and the invariant mass of J {X) is the global momentum and invariant mass of the set 
of jets (invisible particles). In this way, many of the results of the previous subsection 
can be applied here replacing the Vl^-boson and its mass by the two initial supersymmetric 
particles and their global invariant mass, A/inv, or equivalently the total energy at CM; 
and replacing e — t- J, u ^ X. The main difference with the W-case is that J and X 
have masses different from zero, which are equal to the invariant masses of the visible and 
invisible systems respectively. These masses, which change from event to event, may be 
quite large and cannot be neglected. In addition, Minv does not coincide with the sum of 
the masses of the initial supersymmetric particles, since they are not produced exactly at 
rest, although this is normally a good approximation. 

Now, working in the CM system, we can formally consider a subset of events with 
identical structure of jets and invisible particles, differing only in the 0, (p angles at which 
the pseudo-particle J (and thus X) is produced. Obviously for those processes the phase 
space contains a factor like the one given in eq.(2.5), replacing e — )■ J, — )■ X. And this 
holds for any jet configuration. So we expect a similar pole in the histogram. This can 
also be seen by performing a change of variables p ^ ^ p = p ^ in the phase space 
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integration of the jet momenta. One gets 

TT d^P^ d^P^ -TT d^P^ .2im 

(27r)3 2Ei (27r)3 2E^ JLi (2vr)3 2Ei ^ ' ' 

and a similar transformation for the momenta of the invisible particles. Then, in analogy 
with the VF— case, we can now consider the differential cross section with respect to the 
transverse component of the J— momentum, pj., obtaining that eqs.(2.4, 2.5) hold with the 
e — )• J, — X replacements. Thus 



dcos 9 



da oc > (2-11) 



where 

d cos 6 —pj. Ip"^ I 



dPT J\pJ\^ - {pJ^Y 



-2E p^lir'l 

^^'^ ' (2.12) 



_ (^J + E^Y^E'^ - {E^ - E^f 
with E = E"^ + E^ denoting the total energy in the CM and 



= ^{P^r + mj, = ^{E^Y - {pfY ■ (2.13) 

Again, we define a global "transverse energy" as 

£t = E^ + E^, (2.14) 

so that the differential cross section in the transverse momentum, eqs.(2.11, 2.12), shows 
a pole at 

^rlpoic = E = A/i„v ~ mi + 7712 . (2.15) 

Here mi,m2 are the masses of the supersymmetric particles initially created, which, in the 
last identity, we have assumed to be produced approximately at rest in the CM system. 
This pole is maintained when the cross section is displayed in the £t variable because the 
change of variables pj. — t- £t does not introduce any singular behavior. Note that the 
position of the pole is always the one given by eq.(2.15) for any subset of events (with 
any structure of jets and invisible particles) considered. So a global histogram in the £t 
variable must show a peak aX £t = E mi + m2. This is illustrated in Fig. 3, which 
shows the distribution of £t at CM compared with the invariant mass, Minv, of the the 
supersymmetric particles initially produced. The plot corresponds to the benchmark point 
SU9 (defined and discussed in sect. 5) where for this particular example we have not 
included Initial State Shower. As one can see, the two histograms are pretty similar, 
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Figure 3: £t at CM frame (pink) compared with the global invariant mass, Mjnv, of the super- 
symmetric particles initially produced (grey) for the benchmark point SU9 defined in sect. 5. 



exhibiting a peak at the expected value. This similarity is mainly because J and X are 
normally rather heavy and therefore mj and nix are the dominant terms in the invariant 
mass, leading to a small contribution of p^, even for events outside the pole. 

We have indeed checked the presence of the peak at the expected value {£t — "ii + 
1112) for a large collection of CMSSM and more general MSSM models, using a PYTHIA 
simulation. This fact will be used in the next section to understand the correlation between 
Msusy and M^s. We postpone the details and discussion of this and other related checks 
to section 5. 

Since £t is invariant under boosts in the z— direction, the previous pole shows up also in 
the LAB system, provided the net transverse momentum of the two initial supersymmetric 
particles is small, which is the usual case. Once more, if desired, the effect of such net 
non- vanishing transverse momentum can be formally extracted by using, instead of £t, the 
transverse mass, Mt, defined as in eq.(2.9), i.e. = £^ — {pj, +p^)^. 



3. Understanding the correlation Msusy — M^s 

In this section we show that the empirical correlation between Meg (at the maximum of the 
histogram of events) and Mgusy is an echo of the fact, discussed in the previous subsection, 
that the cross section has a pole (in practice a maximum) at £t — "ii + m2 ■ 

Since the correlation has to do with the total mass of the two supersymmetric particles 
initially created, a convenient definition of -Msusy is 

b ^ab {ma + nib) , , 

Msusy = ' ^ (3.1) 

where a, h run over all supersymmetric particles and (Jab is the production cross section of 
the {a, 6} pair. Roughly speaking MgusY is the sum of the masses of the two supersym- 
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metric particles in the dominant channel of production. This definition is exactly twice 
the definition of Mg^sy proposed in ref.[f6] and approximately twice the HPSSY definition 
given in eq.(1.2) if the dominant channel is squark-squark or gluino-gluino. In the rest of 
the paper Ms^sy will always refer to the definition of eq.(3.1)^. 

Let us now discuss the definition and meaning of Meg. Certainly Meg has some 
resemblance with the £t variable defined in eqs.(2.14, 2.13). In order to deepen in the 
connection between both, let us re-write explicit expressions for the two variables. Mefj is 
defined as 

Meff = W\ + \PT^n = Af/g + \pt^'^^\ , (3.2) 
j 

where for the sake of the discussion we have separated the jet contribution to the effective 
mass, denoted by M^^, from the invisible one, Ip^^"^*^!. Normally is by far the most 
important contribution to Meff. Note that we have slightly modified the initial definition of 
eq.(l.l) by extending the sum to all jets, rather than just the four hardest jets. This seems 
reasonable and follows the strategy of ref.[16]. A different thing is how to impose cuts in 
order to optimally count all the jets coming from the partonic process, leaving outside all 
the jets coming from initial state radiation. This will be addressed in sect. 5. 

On the other hand, £t, defined in eqs.(2.14, 2.13), can be re-written in the following 

way 

£t — + , (3.3) 

with 

{E^ f = {E'f - {Pi? 

N N 

= Yl (l^T 1^ + "^j) + Yl V l^rP + ^1\/\Pt\^ + ^1 cosh (vi - yj), (3.4) 
(E^f = '{E^f - {P^f 

2 

= Y (IPt f + "^L) + ^V^Ip/' P + "^li V'IPt + cosh (2/xi - yx2)> (3.5) 

where yj and ^ are the rapidities of the jets and the two neutralinos (or whatever LSPs) 
respectively. 

Comparing the above expressions (3.2)~(3.5) for M^q and £t we see that M^^ < Ej. 
and Ip™*^*^! < E^ , so Mcfj < Et- Focusing in M^^, which is the main contribution to the 
effective mass, we see that the equality -M^ = Ej. is only achieved when the jet masses 
are negligible, which is often a good approximation, and all jets have the same rapidity, 
yi = yj, which is trivially satisfied when there is only one jet, but for two or more jets is 

^All these definitions may become problematic when the the masses of the more frequently produced 
pairs of superparticles are very different. This instance will be discussed in sect. 5. 
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never satisfied in practice. Hence, using instead of Ej, makes the peak of the histogram 
to occur systematicahy below mi + m2- 

We can go further by estimating how large is this effect in a typical case. Working in 
the CM system, the events corresponding to the peak of the histogram have — 0- They 
satisfy Ej, = E^ = Y^- E^ ~ Ylij IP'^ I; where in the last equality we have neglected the jet 
masses. Consequently, at the peak of the histogram 



= ^ . (3.6) 

Usually SUSY searches are done by considering events with a rather large number of jets, 
see e.g. refs.[13, 14, 17, 18]. Hence, for the sake of the estimate, a reasonable simplification 
is that (for the events at the peak of the histogram and working in CM) the directions of the 
various jets are distributed in a more-or-less random way. Then the differential probability 
that a jet occurs at a particular 6 is sinOdO. In average {p>^) = |p J | / d9sm^e= (7r/4)|p J |. 
Therefore, at the peak of the histogram 

l^. = .78.5%. (3.7) 

This correlation is the main reason for the correlation found between Mgfj at the peak of 
the histogram and ~ 80% Msusy 

On the other hand, the contribution to the effective mass in eq.(3.2) is also 

systematically smaller than E^ , defined in eq.(3.5). As we will see in the next section 
and Appendix A, for events where the invisible particles are just two neutralinos with 
momenta larger than their masses, a more fair estimate for E^ is 2|p™''^'^| . This means that 
the invisible contribution to the effective mass is around half the value suitable to get the 
peak of the histogram at mi + m2. 

We have numerically checked that the peak of the Mcfj histogram is indeed around 
70%— 80% Msusy; the precise value depends on the model and the cuts used in the analysis. 
This is illustrated by the statistical survey of CMSSM models presented in Fig. 5 below 
(light green crosses), which will be discussed in more detail in sect. 4. 

4. Proposal of a new kinematic variable, 

The discussion of the previous section not only allows to understand the correlation between 
Mcff and Msusy; it also allows to propose an alternative variable which shows an even more 
robust correlation. As argued above, an histogram in St shows a maximum near Msusy! 
so the idea is simply to devise a new variable which is both measurable and as close as 
possible to St- 

Examining eqs.(3.3, 3.4, 3.5), we see that Ej,, which is the dominant contribution to 
£t, can in principle be extracted directly from experiment. E^ cannot be deduced from 
the experiment, but is clearly larger than Ip™^*^*^! (the quantity used in the definition of 
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Meff), unless the masses of the neutraUnos are fairly smaller than their momenta and the 
latter have exactly aligned directions, which is unlikely. A much better estimate can be 
obtained by assuming that the relative directions of the two neutralinos are random in the 
CM system. This is exact if the the two initial supersymmetric particles are created at 
rest in CM. Under the further assumption that the 3-momenta of the two neutralinos are 
similar in magnitude and larger than the neutralino masses, it turns out that, in average, 
(E^) ~ 21^™'^*^!; for more details see Appendix A. Thus our new kinematic variable, say 
simply reads 



where Ej, is given by eq.(3.4) or eq.(2.13). In the next section we will test the "performance" 
of £^ as a tool to determine Msusy, and compare it to that of Mgff, 

To finish this section, let us mention some of the a-priori advantages (and one disad- 
vantage) of £^ with respect to M^s- The most obvious advantage is that is much closer 
to £t, and so we expect the peak of the corresponding histogram to be quite close to Msusy 
(which is in fact the case, as we will see). Notice that the definition of £lj^ (in particular 
the jet contribution, Ej.) contains information not only about the transverse components 
of the jet-momenta but also about the longitudinal ones, thus being more informative than 
Meff. Actually, M^s is proportional to £t only as an average, and the precise proportional- 
ity factor changes with analysis-dependent features, such as the number of jets considered 
for the selected events or the cuts performed on the various kinematical variables. Con- 
sequently the peak on the histogram is likely to indicate more faithfully the value of 
Msusy than the Meg— one. Another important advantage of £^ is that it is robust under 
the procedure followed to identify the jets. Actually, as it is clear from eq. (2.13) or (3.4), 
for the definition of Ej, we do not even need to talk about jets: we could perform the 
sum directly over hadronic final states; or we could consider all the hadronic particles as 
forming a single jet. The result is the same. On the other hand, a disadvantage of £^ 
with respect to Mgff is that it relies on a good knowledge of the longitudinal components 
of the jet momenta. E.g. one might dismiss a jet which really arises from the partonic 
event because it has large longitudinal component and does not pass the cut in rapidity. 
This would distort the estimated value of This happens also for Mgff, but in the latter 
case the contribution of such jets is less important, since only the transverse component is 
counted. 

In consequence, £!j^ and M^s are complementary variables, rather than competitors; 
and plotting histograms in both can be useful to cross-check the results, allowing a better 
and more robust identification of Mgusy The extension of this procedure to other scenarios 
of new physics (not necessarily SUSY) is also straightforward. 

5. Testing the efficiency of M^s and £^ in assorted SUSY models 

We will discuss now the comparative behaviour of M^g and £^, and other kinematic 
variables, in the context of different MSSM models. We will simulate LHC signals at 14 TeV 



£i 



■eff _ 
T — 



E^ + 2\pt 



miss 




- 11 - 



center-of-mass energy using SOFTSUSY [19] and PYTHIA version 6.419 [20]. Focusing on 
events with multijets + missing transverse momentum, and applying the foUowing cuts 

1. At least three jets with pT > 50 GeV. 

2. The hardest jet with px > 100 GeV and \rj\ < 1.7. 

3. p^'""' > 100 GeV. 

4. A0(jeti -p^''''') > 0.2, A(/.(jet2 -p??^"") >0.2, A(/)(jet3 - p™^^) >0.2. 

For the evaluation of the various kinematic variables we will count all the jets (with pT > 50 
GeV, as mentioned above)^. For the construction of the jets we will use FAST JET [21], 
with the antikt algorithm, E scheme and R = 0.4. 

Let us start by illustrating the behavior of the various kinematic variables by consid- 
ering a typical SUSY model, namely the benchmark point SU9, defined in ref. [12] and 
specified by the following values of the supersymmetric parameters: 

m = 300 GeV, M1/2 = 425 GeV, ^ = 20, tan/3 = 20, ^ > 0. (5.1) 

where m, and A are the universal scalar mass, gaugino mass and trilinear scalar 

coupling (all quantities defined at the Mx scale) respectively; tan/3 = V2/V1 is the ratio 
between the vev's of the two Higgs doublets, Hi, H2; and fi is the usual Higgs mass term 
in the superpotential, defined also at Mx- The corresponding values of the squark mass 
(first two generations) and the gluino mass are rriq = 920 GeV, Mg = 994 GeV. The 
corresponding value of Mgusy, defined according to eq.(3.1), is 

Msusy = 1898 GeV . (5.2) 

Figure 4 (top- left panel) shows the histogram distributions of the invariant mass, Mjnv = 
Ec M , and the histogram distribution of M^s ■ Clearly the Mjnv distribution shows a sharp 
starting point, almost coinciding with the maximum, at Mjnv — -^SUSY- Although the 
peak is quite sharp, the width of the distribution indicates to what extent the initial 
supersymmetric particles are created with non-vanishing momenta at CM. Of course Minv 
is not directly measurable, so this histogram cannot be built in practice. On the other 
hand, the Meff— histogram shows indeed a maximum correlated with MgusY, namely Mcs — 
70% Msusy in this case. Figure 4 (top-right panel) shows the distribution of the £t variable 
compared with the Minv— one. As expected from the discussion of sect. 2, the maximum of 
the two histograms are indeed very close. Comparing this with the histogram, figure 
3, one can see that £t is shifted to the left with respect to £q<^. Note that now we are 
doing a more realistic simulation, including Initial State Shower (ISS), so, in order to avoid 
contamination from these ISS, one has to require jets with high pj- (as described above). 

^We have checked (using a few SUSY points) the stabihty of £!j? when applying a cut in 7;. Once one 
requires pr > 50 GeV, the impact of such cut is very small. The histograms are very similar whether one 
uses a I77I < 3 cut or not cut at all. 
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Figure 4: Minv and Mce (top- left panel), Minv and St (top- right panel), Afinv and (bottom-left 
panel) and 2|pji"'*'''| and (bottom-right panel) for the SU9 benchmark point. 



Hence, the main reason for this shift is that by including those cuts we are loosing some of 
the final particles coming from supersymmetric decays, which leads to an underestimate of 
the visible part; but still the maximum is where was expected. Fig 4 (bottom-left panel) is 
a similar plot, but using the measurable variable, defined in eq.(4.1). We see that the 
£:|?'-histo gram maintains the peak close to Mgusy, though is somewhat less sharp than in 
the 1?^— histogram. This is due to the fact that the invisible contribution to £"1? in eq.(4.1), 
i.e. 2\pr^^^^\, is an average of , i.e. the actual contribution, to £t] see eqs. (3.3, 3.5) 
and Appendix A. The goodness of such average (quite satisfactory in this case) can be 
appreciated in Figure 4 (bottom-right panel), where these two quantities are plotted. 

Let us now explore how good is the behavior of M^g and in general SUSY models. 
We start considering a sample of 500 random CMSSM points requiring dominant qq, gq, gg 
production, which is the usual case. Recall that the CMSSM is defined by the values of 

m, Ml/2, tan/3, sign/x . 

(Incidentally, the benchmark point SU9, defined at eq.(5.1) was a particular CMSSM 
model.) For all the models we have performed the simulation of the proton proton collisions 
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Figure 5: Mgusy versus £^ (blue stars) and Afcff (light green crosses) for the CMSSM. 



using the specifications presented at the beginning of this subsection. 

Figure 5 shows the values of £^ (blue stars) and and M^s (light green crosses), which 
maximize their corresponding histograms, versus Mgusy for those 500 CMSSM models. 

As expected from the previous discussion, there is an remarkable correlation between 
£f at the maximum of the histogram and Msusy] namely £^ ~ Afsusy The correlation is 
also good for the effective mass: Mcfj (at the histogram maximum)~ 70% Msusy, confirming 
the empirical observations of previous literature [15, 16]. 

Notice that, for large values of Mgusy, there are few points which do not correlate well. 
This is because at large energies the production of charginos and neutralinos may become 
competitive. In that case there appear two separate peaks in the histograms, which makes 
the definition of Mgusy fairly contrived. 

Let us now extend the previous study to more general MSSM models, allowing non- 
universal soft parameters at the Mx scale. For that goal, we extend the previous parameter 
space, eq.(5), to a 15-dimensional parameter space, defined by 

Mi,M2,M3; At,Ab,Ar; mH^^H/, mq^^g^,m^^ -^^,m^^,mi^; m^-^i^, "if^.i/.^ , m^^; tan/3. 

Here Ma are the (non-universal) gaugino masses; and Ai, rrii are the non- universal trilinear 
scalar couplings and scalar masses {i denotes flavour species), I denote squarks and 
sleptons of the first two generations. As for the CMSSM case, we have studied 500 models 
chosen at random in this parameter space. 

Figure 6 (left panel) is as Figure 5 but within this more general MSSM scenario. In 
this case the correlation is not as good as for the CMSSM, something that was empirically 
noticed in ref.[16] for the Mgg variable. We have checked that the reason is that now 
the neutralino masses can be much larger than in typical CMSSM models. Then, since 
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Figure 6: Left panel shows Mgnsy versus £^ (blue stars) and M^s (light green crosses). Right 
panel shows Msusy versus £t (red crosses). Both for the MSSM. 



in the definition of both Mefj and f|? one is neglecting the neutralino masses, one is 
missing a potentially important piece. Furthermore, in some cases many neutrinos can be 
produced. This makes the average missing transverse energy, ^ larger than the average 
~ 2|p™'^''| used in the definition of £^ (and of course larger than the piece entering 

the definition of Mgff). This can be checked by plotting the true £t variable (which cannot 
be directly measured). This has been done in the right panel of figure 6. Remarkably, here 
the correlation with Msusy is nicely maintained, indicating that the spreading exhibited by 
and Mcff is indeed due to an underestimate of the invisible contribution. 

The previous discussion suggests that a better correlation could be found by refining 
the estimate {E^) ~ 2|p™''^'^| used above, by taking into account the finite size of neutralino 
masses. As discussed at the end of the Appendix, a more refined estimate of the invisible 
transverse energy for the events at the pole is obtained by retaining the dominant terms in 
m^, and reads {E^) ~ ^4:{p^)'^ + 4m2 ~ 2p^ + 2m'^/ E^ . Using the fact that typically 
E^ contributes less than half to £t, we get a conservative correction ~ 4m^/Msusy to £^ 
for the events at the pole. Alternative, we can keep the definition of £^ given at eq. (4.1). 
Then the value of £^ at the maximum of the histogram must be around 



M 



susy 



M 



susy 



(5.3) 



rather than Msusy Incidentally, a similar correction to Msusy with the same functional 
dependence but half the value, was proposed in ref.[16]. 

2 

Fig. 7 shows Msusy— 4 j^^^^^ versus £f for the same 500 general MSSM models displayed 
in Fig. 6. Clearly the correlation is now much better, comparable to that found for CMSSM 
models (Fig. 6). The correction (5.3) can also be applied to CMSSM models, but in that 
case the improvement is less significant since is normally much smaller than Msusy 
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Figure 7: Afgusy — 4 ^^^ ^ versus £|? (blue stars) and M^s (light green crosses) for the MSSM. 
6. Conclusions 

There is an intense activity to explore assorted strategies for the search of new physics 
at the LHC. A paradigmatic case is supersymmetry (SUSY). Under the assumption that 
superparticles are light enough to be produced and detected at LHC, one can expect a 
typical signal of multijets and large missing energy. However, the analysis is quite compli- 
cated and the correct identification and reconstruction of the properties of the final states 
is quite tough. Besides, the separation of the decay chains is normally not possible. In this 
context the use of appropriate kinematic variables is instrumental to maximize the signal 
of new physics from the Standard Model background. 

A handy and extensively used kinematic variable, e.g. in ATLAS analyses, is the 
Effective Mass variable, defined as 

Meff = ^ \pS\ + \PT^n = M/ff + \pt^'^^\ . (6.1) 
j 

It was empirically found in refs.[15, 16] that for the Constrained Minimal Supersymmetric 
Standard Model (CMSSM) there is an outstanding correlation between the value of M^s 
at which the histogram of events has a maximum and the typical supersymmetric masses. 
Namely, 

^effUax ^ 80% Msusy , (6.2) 

where M^usy ~ mi+m2, i.e. the sum of the masses of the supersymmetric particles initially 
created (usually squark-squark, gluino-gluino or squark-gluino) . However, for some points 
in the CMSSM, and especially in the general MSSM, the correlation fails. 

In this paper we have found an explanation for the above correlation (6.2). We have 
argued that Meff is typically 80% of a kinematic variable, £t, defined as the sum of the 
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transverse energies of the visible and the invisible parts of the decay products of the initial 
supersymmetric particles; for more details see eqs.(2.13) and (3.4, 3.5). Then it can be 
shown that, at CM, the cross section has a maximum at £t = -E'cM ~ fni + ^2- This 
relation is invariant under longitudinal boosts and thus is kept as a good approximation in 
the LAB frame. 

This understanding has allowed us to propose a new kinematic variable, the effective 
transverse energy 8^ ^ which is measurable and as close as possible to £t- It reads 

£f = E:^ + 2\prlP-^^\ (6.3) 

where (Ej,)"^ = [E'^)^ — (p^)^ is the hadronic transverse energy [for more details see 
eqs.(2.13), (3.4)]. £^ is a kinematic variable alternative to Meg, which shows an even 
better and more direct correlation with the supersymmetric mass 



£: 



eff 



Msusy . (6.4) 



Besides, £^ has other advantages, like being more robust under the procedure followed 
to identify the jets (actually, it is completely independent of it). On the other side, the 
determination of £^ relies (more than for M^g) on a good knowledge of the longitudinal 
components of the jet momenta. In consequence, £^ and Mcs are complementary vari- 
ables, rather than competitors; and plotting histograms in both can be useful to cross-check 
the results, allowing a better and more robust identification of Mgusy The extension of this 
procedure to other scenarios of new physics (not necessarily SUSY) is also straightforward. 

We have shown these features in the context of the CMSSM or more general MSSM 
models by examining the above correlations in a large number of points in the parameter 
space. The results show that £^ is a simple and valuable kinematic variable to determine 
the characteristic supersymmetric masses. 

Our analysis also shows why the above correlations fail for some models. This is 
because at large energies the production of charginos and neutralinos may become compet- 
itive. In those cases there appear two separate peaks in the histograms, which makes the 
previous analysis too simple. In addition, for general MSSM models (departing from the 
conventional CMSSM) the neutralino masses can be much larger than in typical CMSSM 
models. Furthermore, in some cases many neutrinos can be produced. In those cases, the 
term representing the missing transverse energy in (6.1) and even in (6.3) underestimates 
the actual invisible contribution and a more refined estimate of it must be used. Then, re- 
taining the dominant contributions in the neutralino mass, the right-hand side of eq. (6.4) 
gets shifted as Afsusy Msusy — 4m^/Msusy, significantly improving the correlation. 

7. Appendix 

In this appendix we evaluate the average value of the missing transverse energy, defined in 
eq.(2.13), 

{E^? = {p^f + ml = {E^f-{pff. (7.1) 
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Here , pf are the transverse and longitudinal components of the total invisible 4- 
momentum 



X 



P 



P1+P2 , 



(7.2) 



where pi,P2 are the two 4-momenta of the two final-state neutralinos, XiiX2- From now 
on we will work in the CM reference frame, see Fig. 8. As represented in Fig. 8, 9 denotes 
the angle oi p ^ with the longitudinal axis, z; and a the angle between pi and p2. Let us 
now write the missing energy, , using the second equality in eq.(7.1): 



= 2ml + (Pi + ^"^\)iP2 + "^x) + + P-^) '^"^^ ^ ~ 2|pi| |p2| cos a cos^ 6^. (7.3) 
On the other hand, the invisible transverse momentum, p^, is given by 

{p^ f = {p^ + p^) sin^ e + 2\pi I |p2| cos a sin^ 6* . (7.4) 



Hence, 



{E^f = {p^f + 2ml + 2^{p^^ + ml){p.2^ + ml)-2\pi\\p2\cosa . (7.5) 

This expression admits further simplifications under some reasonable assumptions. Assum- 
ing m\ <^ {Pi)'^, we have 



{E. 



X\2 



{p^f + 2\pi\\p2\{l-cosa) = (p#) 



X\2 



1 + 



2|pi||p2|(l - COS a) 



sin^ 6 [p-^ + P2 + 2\pi I \p2 \ COS a) 



(7.6) 



Finally, we can assume that typically \pi \ ~ \p2\- Then 



{E^f^ip^ 



X\2 



1 + 



1 — COS a 



sin 0(1 + cos a) 



(7.7) 
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Now, we can single out the previous expression for the events at the pole, i.e. with 9 = 7r/2. 
Notice also that, assuming that the initial supersymmetric particles are approximately at 
rest, it turns out that the relative directions of pi and p2 are random, so the differential 
probability for a particular value of the angle a is V{o) = ^ sin a da. In conclusion 

1 

{E^)po\e - i^Pt I da sin a 
^ JO 

which is the estimate used of throughout most of the paper. The goodness of this 
approximation is illustrated in Fig. 4 (bottom-right panel) for a typical example. 

A more refined estimate can be obtained by keeping the dominant terms in in 
eq. (7.5). Then, eq. (7.8) gets slightly modified: 

{E^ )poic ^ ^4(p^)2 + 4m2 . (7.9) 
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